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I. INTRODUCTION 

In nuclear physics, the pairing mechanism is one of the 
basic ingredients in determining the nuclear structure prop- 
erties 1 1 ], and becomes even more significant for the exotic 
and superheavy nuclei which have attracted wide interests 
of the community during the past decades Ii2 ffl"2l . Specifi- 
cally in exotic nuclei [13], the weakly bound nuclear systems 
under extreme condition, crucial roles played by the pairing 
correlations are found not only in predicting the isospin lim- 
its of finite nuclei but also in developing the exotic 
modes, such as nuclear halo structures and possible BCS-BEC 
crossover 0^9). For superheavy nuclei with extra large charge 
numbers, not only the detailed microscopic structure, but also 
the bulk properties, e.g., nuclear shapes, fission barriers and 
collective modes, are essentially related with the effects of 
pairing 

While experimentally it is not so straightforward to mea- 
sure the pairing effects, which are in general evaluated by 
the odd-even staggering (OES) on the nuclear binding en- 
ergy. The binding energy of a system with odd nucleon (neu- 
tron or proton) number is found to be lower than the arith- 
metic mean of two even neighbors, which leads to the so- 
called OES of single-nucleon separation energies along the 
isotopic or isotonic chains. Since the early days of nuclear 
physics, the OES was interpreted as the presence of pairing 
correlations between nucleons in nucleus (TJ. In the indepen- 
dent quasi-particle picture, the OES extracted from the exper- 
imental binding energies is often taken as the reference of the 
pairing-gap energy, which also provides a quantitative observ- 
able to constrain the pairing interaction Ifl4ll . 



During the past decades, many successes have been 
achieved by the covariant density functional (CDF) theories 
in describing the structure properties of nuclear systems ITT3I — 
1211 . One of the most popular CDF models is the relativis- 
tic Hartree approach with the no-sea approximation, namely, 
the relativistic mean field (RMF) theory If22l l23l . which has 
been widely applied in exploring the properties of both ground 
states G6] [H [HI and excited states [18] |24) for the finite 
nuclei in and far from the valley of [5 stability. However, 
due to the limit of the approach itself, significant system de- 
grees of freedom are missing in RMF, such as the one-pion 
exchange ll25l and p-tensor couplings. With the growth of 
computational facilities and the development of new meth- 
ods, such defects can be eliminated with the inclusion of ex- 
change (Fock) terms, which leads to a new CDF model — the 
density-dependent relativistic Hartree-Fock (DDRHF) theory 
CUl and its natural extension, the relativistic Hartree-Fock- 
Bogoliubov (RHFB) theory [26|. Besides compatible quanti- 
tative accuracy as RMF in describing nuclear bulk properties, 
substantial improvements are also achieved by DDRHF in the 
self-consistent description of nuclear shell structures and the 
evolutions ET1 l27l 1281 . the relativistic symmetry restoration 
C91 - I3T1 , low-energy excitation modes [32], and neutron star 
physics |[33][33 etc. 

In general the pairing effects in open-shell nuclei are con- 
sidered within the BCS or Bogoliubov schemes. For a stable 
nuclear system, the BCS method can provide an efficient and 
simple way to handle the pairing correlations while it meets 
serious problems going beyond the stable region, especially 
when the halos emerge fl9l . Approaching the nuclear isospin 
limits, the single neutron or proton separation energy becomes 
comparable to the pairing-gap energy such that the continuum 
can be got involved easily by the pairing correlations. In terms 
of Bogoiubov quasi-particles, both mean field and pairing cor- 
relations are unified into the Bogoliubov scheme such that the 
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continuum effects can be naturally taken into account. Aiming 
at the weakly bound nuclear systems, the CDF models com- 
bined with Bogoliubov framework have been extended as the 
relativistic Hartree-Bogoliubov (RHB) theory |[T8l[35l and the 
RHFB theory l26l . Compared to the former, the RHFB the- 
ory provides a more concrete platform to evaluate the pairing 
effects due to the improvements brought by the Fock terms, 
especially with the inclusion of p-tensor couplings BTTl . 

Besides the self-consistent treatment of the continuum ef- 
fects, e.g., by the Bogoliubov theory, the theoretical reliability 
in describing the pairing effects also depends on the adopted 
pairing interaction. In the realistic calculations the particle- 
particle (p-p) pairing interactions are usually taken as a phe- 
nomenological form, such as zero-range ^-forces or finite- 
range Gogny interactions with great success in the relativistic 
and non-relativistic calculations [19, 26l l36l [37ll . Due to its 
numerical simplicity, the <S-type forces are popularly applied 
in dealing with the pairing effects, especially when a realistic 
density dependence is introduced 08l . While limited by the 
zero-range formalism, it may require a sophisticated energy 
cut-off to decide the pairing window, as well as the strength 
of pairing interaction [38]. In contrast to the simple (5-force, 
the finite-range Gogny forces have been widely taken as the 
effective interactions not only in the p-p channels but also in 
the particle-hole (p-h) channels and potentially better system- 
atics is expected for the p-p interaction. In addition benefited 
from the characteristic finite range, a natural energy cut-off 
is already embedded in modeling the pairing space. While it 
should be noticed that the pairing parts of Gogny forces were 
adjusted in consistent with the non-relativistic mean fields 
071 [39l l40l . The global consistence associated with the rela- 
tivistic scheme BT1I421 is still required to be tested, especially 
when the Fock terms are included in the p-h channels. 

In this work, the OESs on nuclear binding energy will be 
systematically studied under the CDF scheme, specifically the 
RHFB and RHB theories. The content is organized as follows. 
In Sec. [Il] we briefly introduce the general formalism of the 
RHFB theory with finite-range Gogny pairing force. The neu- 
tron OESs along the isotopic chains of C, O, Ca, Ni, Zr, Sn, 
Ce, Gd and Pb, and the proton ones along the isotonic chains 
of N = 50, 82 will be discussed systematically in Sec. Ill 



as 



well as the effects with modified pairing strength. Finally, the 



summary is given in Sec. IV 



II. 



THEORETICAL FORMALISM AND NUMERICAL 
DETAILS 



Relativistically the nucleon-nucleon interaction in nuclear 
system is mediated by the exchange of mesons and photons. 
Standing on this criterion, the CDF model Lagrangian, i.e., the 
theoretical starting point, is constructed by including the de- 
grees of freedom associated with the nucleon (if/), the isoscalar 
cr- and w-mesons, the isovector p- and 7r-mesons, and the pho- 
tons (A) 11431 . Following the standard variational procedure, 
the system Hamiltonian H is then determined as well as the 
field equations for nucleons, mesons and photons, respectively 
the Dirac, Klein-Gordon and Proca equations Il23ll43l . 



Standing on the level of mean field approach, the contri- 
butions of the negative energy states are generally neglected, 
namely the no-sea approximation. The CDF energy functional 
is then determined by the expectation of the system Hamilto- 
nian H with respect to the Hartree-Fock ground |i>o), 



E =<0 |//|cDo>, |O > =Y[ c] |0) 



(1) 



where the index i denotes the states with positive energy and 
|0) is the vacuum state. Within the RHB theory, the en- 
ergy functional ([T]i contains only the Hartree contributions, 
whereas both Hartree and Fock terms are included explicitly 
in the RHFB theory EH |35]. 

In the spherical system, the variation of the energy func- 
tional ([T]l with respect to the Dirac spinor tf/{r) leads to the 
relativistic Hartree-Fock (RHF) equation as 



dr' h(r,r')i[r(r') -eij/(r), 



(2) 



where e is the single-particle energy and h(r, r') denotes the 
single-particle Hamiltonian J26). Combining the Bogoliubov 
transformation B4l with CDF models, the above RHF equa- 
tion can be extended into the RHFB one |45 1 as 



I 



dr' 



h(r,r') A(r, r') 
-A(r, r') h(r,r') t 
(A + E, 



( iffu(.r') 











A-E, 



Mr) 



(3) 



where \pu and i/<y are the quasi-particle spinors, E q denotes 
the quasi-particle energy, and the chemical potential A is intro- 
duced to preserve the particle number on the average. In the 
single-particle Hamiltonian h(r, r'), the retardation effects are 
neglected as is usually done in mean-field calculations 11261 . 
The pairing potential A(r, r') in Eq. ^ reads as 



A„(r, r') = - - ^ V'^(r, r')Kp(r, r'), 



with the pairing tensor 

K a (r,r') =ilf Va (r)*il/ Ua (r'). 



(4) 



(5) 



For the p-p interaction V pp in Eq. Q, the following finite- 
range Gogny force is adopted in this work, 

V(r, r') = V e ((r - r ' )lp > )2 (Wi + BjP a - HiP T - M,P' r P T ), (6) 

(=1,2 

where , Wj, B h Hj, and Mj are the parameters of Gogny pair- 
ing force. 

Except for few special cases, e.g., the RHB model with 
zero-range pairing force, the RHFB equation Q is generally 
in an integro-differential form, which is difficult to be solved 
efficiently in the coordinate space. In this work, the Dirac 
Woods-Saxon (DWS) basis l46l . which can provide appropri- 
ate asymptotical behavior for the continuum states, is intro- 
duced to solve such integro-differential equation. For further 
details, the readers are referred to Ref. [26|. 
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To evaluate the pairing effects, the following three-point in- 
dicators [47, 48] of the OES of nuclear binding energies are 
introduced for isotopes (v) and isotones (tt), respectively, 

(-D N 

A< 3) (A0 = y —^- [S„(N,Z) - S n (N + 1,Z)] , (7a) 

(Z) [S P (N,Z) - S P (N,Z+ 1)] , (7b) 

where A v (N) and A„ (Z) represent neutron (v) and proton 
(tt) OES, i.e., the difference on single-nucleon separation en- 
ergies S n (or S p ) of two neighboring isotopes (or isotones). 
Quantitatively there also exist other higher order indicators to 
evaluate the pairing effects, such as the four-point formula, or 
the five-point ones 11471 

At 5) (A0 = ( -^[S n (N + 2,Z) - 3S n (N + 1,Z) 

+ 3S n (N,Z)-S n (N- 1,Z)], (8a) 

A«(Z) =t^-[S p (N,Z + 2) - 3S P (N,Z + 1) 

+ 3S P (N,Z)-S P (N,Z- 1)], (8b) 

However, these higher order indicators may get the mean-field 
features involved, which may partially smooth out the odd- 
even oscillations due to the pairing effects 0481 . Even with the 
three-point one the obtained OESs for the isotopes (iso- 
tones) with even neutron (proton) number are still sensitive to 
both pairing and mean-field effects [48-50]. In order to avoid 
the disturbances of the mean field as much as possible, here 
we only take the OESs [see Eq. ([7]i] of the isotopes (isotones) 
with odd neutron (proton) number as the observables to study 
the pairing effects. 

To have an overall understanding on the pairing properties, 
we performed the calculations within both RHFB and RHB 
for the isotopes of C, O, Ca, Ni, Zr, Sn, Ce, Gd and Pb, and 
the isotones of N = 50 and 82, from which are extracted the 
neutron and proton OESs by Eq. (|7]) for the odd isotopes 
and isotones, respectively. In all the calculations, the p-p in- 
teraction V pp is taken as the finite-range Gogny force D1S 
lIBTI . In the mean field channels, we adopt four CDF effective 
interactions, specifically PKA1 [21 1 and PK03 [27 1 within 
RHFB, and PKDD (52) and DD-ME2 [53 1 within RHB. For 
the selected isotopes and isotones, it is checked to be accurate 
enough to take the parameters of the DWS basis as R mRX = 20 
fm, Nf = 28 and No = 12, where R max is the size of the 
spherical box and Nf{Nd) corresponds to the numbers of pos- 
itive(negative) energy states included in the basis expansion. 
For the nuclei with odd neutron and/or proton numbers, the 
blocking configurations are utilized here to approximate the 
odd-nucleon effects. Practically by blocking different orbits 
around the Fermi surface the configuration with the strongest 
binding is determined as the ground state. Notice that in the 
odd nuclei the time reverse symmetry is broken due to the odd 
nucleon. While the time-odd mean fields are poorly known 
and their effects differ from model to model. In this work we 
neglect the current effects since the mean field would not be 
essentially changed by the odd particle l54ll . 



m. RESULTS AND DISCUSSIONS 
A. Odd-even staggering along the isotopic and isotonic chains 

We first study the neutron OESs of the selected isotopes 
with odd neutron number within the RHFB and RHB theo- 
ries. Figure [7] shows the evolution of neutron OESs A®\N) 
(in open symbols) along the isotopic chains of C, O, Ca, Ni, 
Zr, Sn, Ce, Gd and Pb, where the experimental data (in filled 
squares) extracted from Ref. Il55l are presented as the refer- 
ence. The theoretical results are provided by the calculations 
with PKA1, PK03 and PKDD. The results from DD-ME2 are 
omitted because they show similar systematics as PKDD on 
the OESs. It is clearly seen that all the calculations with the 
selected effective interactions can provide appropriate overall 
agreement with the data to certain extent. 

For the light nuclei, PKA1 shows better agreements with 
the data than PK03 and PKDD on the isospin evolution of 
OESs along the isotopic chain of C. For O isotopes, the ex- 
perimental depressions of neutron OESs A^ 3) (A) at N — 7 and 
N = 15 are well reproduced by all the theoretical calcula- 
tions, consistent with the occurrences of the shell structures 
at N = 8 and 16. While all the theoretical calculations sys- 
tematically underestimate the OESs Ay 3) (A) for both C and O 
isotopes, especially at the proton rich sides. 

For the medium-heavy nuclei it is found that the neutron 
OESs of Ca isotopes are precisely reproduced by PKA1, 
PK03 and PKDD only except the one at N = 31. For Ni 
isotopes, PKA1, PK03 and PKDD also present proper overall 
agreements, except N - 29 for PK03 and PKDD, and N = 39 
for all the selected effective interactions. From the consis- 
tent relation between the pairing effects and shell structures, 
it seems that the shell effect around N - 28 is underestimated 
by PK03 and PKDD, whereas the one at N — 40 is under- 
estimated by PKA1 and overestimated by PK03 and PKDD. 
For Zr isotopes, the neutron OESs are described precisely by 
PKA1 until the distinct discrepancies occur beyond N - 59, 
where the nuclei may be deformed ll56ll . In contrast PK03 and 
PKDD present remarkable deviations from the data in a fairly 
wide range since N > 49. The model deviations in reproduc- 
ing the neutron OESs here may be connected with the fact that 
PKA1, in which the p-tensor coupling is included [21 J, pro- 
vides better descriptions on the nuclear shell structures and 
therefore improves the consistence between the pairing and 
mean-field effects, as mentioned above. 

For heavier nuclei, e.g., the Sn isotopes, the experimental 
OESs At 3) (A) keep around the value of 1.2 MeV, except at 
N - 65 and N = 83. Quantitatively PKDD provides the best 
agreement with the data while the neutron OESs are slightly 
overestimated by PKA1 and PK03. For Ce and Gd isotopes, 
although many of them may be deformed [56 1, the spheri- 
cal calculations with PKA1, PK03 and PKDD still provide 
appropriate descriptions on the neutron OESs, except for Gd 
isotopes beyond N = 89. Along the isotopic chain of Pb, the 
systematics of neutron OESs are properly described by PKA1, 
PK03 and PKDD, while there still exists some systematical 
overestimation. 
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FIG. 1: (Color online) Neutron OESs A|, 3, (A0 [see Eq. (7}] as functions of neutron number N for the C, O, Ca, Ni, Zr, Sn, Ce, Gd, and Pb 
isotopes. The results are extracted from the calculations of RHFB with PKA1 (in open circles) and PK03 (in open up-triangles), and RHB 
with PKDD (in open down-triangles), in comparison with the data [55 1 (in filled squares). The finite-range Gogny force D1S is utilized as the 
effective pairing force. See the text for details. 



Concerning the evolutions of proton OESs along the iso- 
tonic chains of N = 50 and 82, similar systematics are also 
obtained from the CDF calculations with PKA1, PK03 and 
PKDD, as referred to the data [55 1. From Fig. [2] (a) and (c) 
it is seen that the proton OESs for N - 50 and 82 isotones 
are reproduced in an appropriate agreement with data by the 
selected effective interactions. In particular, the calculations 
with both PKA1 and PK03 present the clear depression of 
A®(Z) at Z = 39 for N = 50 isotones, corresponding with 
the sub-shell structure Z = 40, which may imply the extra 
role of Fock terms especially the embedded tensor effects in 
determining the nuclear shell structure for the selected iso- 
tonic chain [57|. Compared to the data quantitatively, how- 
ever, the pairing effects are also somewhat overestimated for 
the selected isotones. It is worthwhile to mention that the cal- 
culations with other relativistic effective Lagrangian, e.g., the 
RHB with DD-ME2 [53], also show similar trend as above. 
In brief one may reach the point that the OESs, coherently the 
pairing effects, are somewhat underestimated by the CDF cal- 
culations for light nuclei and overestimated for heavy nuclei. 

In general the statistical observables, i.e., the root-mean- 



square (rms) deviation o" lms and the average one D from the 
data, are introduced to test quantitative precision of the the- 
oretical calculations. Table [I] shows the values of cr rms and 
T> for the selected isotopes and isotones. It is found that the 
calculations with PKA1, PK03 and PKDD present negative 
average error T> for C and O isotopes and positive ones for 
nuclei heavier than Ni isotopes. Such systematics may also 
indicate that the pairing effects are somehow underestimated 
for light nuclei and overestimated for heavy nuclei by taking 
the finite-range Gogny force D1S as the effective pairing inter- 
action, which has already been demonstrated from the results 
in Fig. [T]and Fig. |2](a, c). 

To eliminate such systematical discrepancy with the data, 
here we consider a simple modification by introducing an ad- 
ditional strength factor £ into the pairing interaction, similar 
as in Ref. BTH , 

V™(r,r') = (V™(r,r'), (9) 

where V pp (r, r') and V^ t (r, r') are the original and modified 
pairing force. The strength factor f is optimized with respect 
to the root-mean-square deviation cr rms from the PKA1 calcu- 
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TABLE I: The root-mean- square deviation cr rms and the average error D of OES for the isotopes C, O, Ca, Ni, Zr, Sn, Ce, Gd, Pb, and for 
the isotones N = 50 and 82, determined by the CDF calculations with PKA1, PK03 and PKDD. The results inside/outside the parenthesis are 
presented by the calculations taking the Gogny force D1S with/without pairing strength factor f as the effective pairing interaction. See the 
text for details. 
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0.90 
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FIG. 2: (Color online) Proton OESs A®(Z) of N = 50 and 82 
isotones as functions of proton number Z. The results are calculated 
by the RHFB theory with PKA1 and PK03, and by the RHB theory 
with PKDD, where the Gogny force D1S is utilized in the pairing 
channel. The experimental data extracted from Ref. [55| are also 
shown for comparison. The results without and with the effective 
pairing strength factor f are shown in the upper ((a) and (c)) and 
lower ((b) and (d)) panels respectively. 



lations with different values of £ since PKA1 can better de- 
scribes the nuclear structure properties [21]. As the exam- 
ples, Fig. [3] shows the neutron OESs of Sn and Pb isotopes 
calculated with the optimized Gogny pairing force D1S, in 
which the strength factor ( are respectively taken as 0.95 and 
0.90. Comparing to the results with original pairing force (i.e., 
£ = 1.0) in Fig. [T](f) and (i), one can see that the agreements 
with the data are remarkably improved with the inclusion of 
the strength factor, where the overestimation on the neutron 
OES is reduced remarkably. In Ref. [42 1 similar conclusions 



have also been reached in reproducing the experimental data 
of inertia moments for very heavy nuclei. For the N - 50 and 
82 isotones, as shown in Fig. |2](b) and (d), the systematical 
overestimation on proton OESs can be also eliminated by re- 
ducing the strength factors as ( — 0.95 and 0.90, respectively. 

In fact, such improvements can be demonstrated quantita- 
tively from the statistical observables in Table [I] for all the 
selected isotopes and isotones. It is shown that with the in- 
clusion of strength factor almost all the CDF calculations 
present nearly zero average deviations D, especially for the 
heavy isotopes and isotones, which corresponds to a system- 
atic improvement on the description of the pairing effects. As 
shown in the second column of Table [I] the strength factors 
£ represent distinct Z-dependence for isotopes, as well as N- 
dependence for isotones, i.e., £ decreases monotonically with 
increasing proton number Z of the isotopes or neutron number 
N of the isotones, consistent with the systematic deviations as 
seen from Fig. [T] and [2] As seen from the rms deviations in 
Table [I] remarkable improvements on the agreements can be 
also found with the inclusion of the strength factor in p-p in- 
teraction. For the neutron OESs the statistic qualities of the 
agreements, i.e., the values of cr lms , are improved by about 
10% for light nuclei (C and O) and about 30% for heavy ones 
(Sn, Ce, Gd and Pb). Whereas for the isotones of N = 50 
and N = 80 (the last two rows in Table [j}, the RHFB calcula- 
tions with optimized pairing forces respectively present about 
30% and 50% improvement on the rms deviations cr rms , which 
keep almost unchanged in the RHB calculations with PKDD 
as well as with DD-ME2. 

As an implement, Fig. [4]presents the proton effective pair- 
ing gaps, i.e., the average gap energies extracted from the 
RHFB and RHB calculations, for the N - 50 (left plots) and 
82 (right plots) isotones with even proton numbers. For com- 
parison, are also shown the OESs of the experimental binding 
energies by three-point (filled squares) and five-point (filled 
circles) indicators. As the experimental reference of pairing 
effects, one may find distinct deviations between these two 
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FIG. 3: (Color online) The same as Fig. [TJbut for Sn (panel (a)) and 
Pb (panel(b)) isotopes calculated by the optimized Gogny pairing 
force D1S with the effective pairing factor f. See the text for details. 
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FIG. 4: (Color online) The average proton pairing gap of N = 50 
and N = 82 even isotones compared with the corresponding neutron 
OES extracted from the data [55 1. The results without and with the 
effective pairing strength factor £ are shown in the left ((a) and (c)) 
and right ((b) and (d)) panels respectively. See the text for details. 



indicators along the isotonic chain of N = 50 [see Fig. [4] 

(a) ] and such deviations become smaller for N = 82 isotones. 
This is mainly due to the fact that these evaluations of pairing 
effects for even-even nuclei have already got the mean-field 
effects involved to different extents [48 49 1. As referred to 
the experimental values, the calculated effective pairing gaps 
still show appropriate agreements, especially for the calcula- 
tions with optimized pairing force. It should be noticed that 
along the isotonic chain of N = 50 the calculations with PKA1 
present distinct depression at Z = 40, which implies the oc- 
currence of some sub-shell closures, also demonstrated by the 
OES depression of the odd isotones at Z = 39 [see Fig. [2] (a) 
and (b)]. While the calculations with PK03 and PKDD do 
not show similar consistence between the results of odd and 
even isotones. E.g., for the odd N = 50 isotones the OES re- 
sults show certain depression at Z = 39 [see Fig. [2] (a) and 

(b) ] whereas the effective pairing gaps for the even isotones 
are just changed smoothly at Z = 40 [see Fig. ffl(a) and (b)]. 
Along the isotonic chain of N - 82, the inconsistence are also 
observed at Z = 58 from the results of PK03 and PKDD, i.e., 
the values of the effective pairing gaps clearly indicate the ex- 



istence of the artificial shell closure Z = 58 ETl 1581 while 
such spurious shell is not indicated by the systematics on the 
OESs of the odd isotones. From this aspect one may find that 
there exists some inconsistence between the relativistic effec- 
tive Lagrangians (PK03 and PKDD) and non-relativistic pair- 
ing treatment. 



B. Nuclear structure properties with optimized pairing force 

Since substantial improvements in reproducing the neu- 
tron/proton OESs have been achieved with the optimized pair- 
ing force y^p t [see Eq. d91], it is worthwhile to test the rele- 
vant effects in describing tne nuclear structure properties, such 
as the single-particle configurations, pairing energies and nu- 
clear binding energies. As a direct response, the occupation 
probabilities of valence orbits will be essentially changed with 
the modifications of the pairing strength, as well as for the 
pairing energies. 

In Table III] are shown the calculated occupation probabil- 



ities of neutron valence orbits 
lhu/2 in 112 Sn, 114 Sn. 
data taken from Refs. 



2^5/2.1^7/2,3*1/2,2^3/2 and 
1I8 Sn and 124 Sn, as compared to the 
The bracketed numbers are 



the results calculated with the optimized pairing force V?? t 
(£ = 0.95). Among the selected effective Lagrangians, PKA1 
shows the best agreements with the data within the error bars 
(denoted by bold types). Whereas less satisfied agreements 
are provided by the calculations with PK03 and PKDD, espe- 
cially the improper order of the pseudo-spin partners 2ds/2 and 
lgi/2- As seem from the data of occupation probabilities, the 
neutron orbit 2d^/2 is almost fully filled whereas lgj/2 is just 
gradually occupied, which implies that 2d$ /2 is bound deeper 
than lgj/2 in the selected isotopes. Evidently PKA1 presents 
more reliable description on nuclear structures than the other 
selected effective Lagrangians. 

In the PKA1 calculations the optimized pairing force V^L 
also bring some systematic improvements on the agreement 
with the data, i.e., being closer to the central values of the data. 
Whereas with PK03 and PKDD, such systematical improve- 
ments can not be observed, which may be partially due to the 
inappropriate order of the valence orbits. If comparing the 
occupations of another pseudo-spin doublet 3 s 1/2 and 2dj/2, 
one may also find the improper order of these two states de- 
scribed by PK03 and PKDD in 112 > 118 - 124 Sn, and by DD-ME2 
in 112 Sn. Compared with the distinct improvements on the 
OESs with VfL the corresponding effects on the occupations 
are relatively weak because the configurations are determined 
not only by the pairing effects, but also more essentially by 
the concrete shell structure, in which the mean field plays the 
dominant role. 

As another direct effect, the modification on the pairing in- 
teraction will substantially change the pairing contributions 
to the energy functional. Taking the even Pb isotopes as the 



representatives, Table III shows the pairing energies from the 
RHFB and RHB calculations with the selected effective pa- 
rameters. For comparison the results calculated by RHB with 
NL1 and by non-relativistic HFB with Gogny D1S are also 
shown [41 1 . It should be mentioned that in the former RHB- 
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TABLE II: Occupation numbers of valence neutron orbits in 1 12,1 14,1 i8,i24 Sn ^ ca i cu i ate( j by RHFB with PKA1 and PK03, and by RHB with 
PKDD and DD-ME2, in comparison with the data |59 60 1 . In the calculations, the pairing force is adopted as the original Gogny D1S and the 
optimized one with ( = 0.95 (inside the parenthesis), respectively. The bold types denote the theoretical results within the error bars of the 
data. See the text for details. 



Orbit 








U2 Sn 












U4 Sn 






Exp. 


PKA1 


PK03 


PKDD 


DD-ME2 


Exp. 


PKA1 


PK03 


PKDD 


DD-ME2 


v2d 5/2 


0.93 


±0.12 


0.82(0.84) 


0.61(0.63) 


0.61(0.62) 


0.65(0.66) 


0.97 


±0.09 


0.86(0.88) 


0.73(0.76) 


0.759(0.80) 


0.80(0.83) 


vlg 1/2 


0.63 


±0.13 


0.55(0.56) 


0.83(0.86) 


0.89(0.91) 


0.86(0.88) 


0.69 


±0.15 


0.65(0.67) 


0.88(0.90) 


0.923(0.94) 


0.91(0.94) 


v3s l/2 


0.24 


±0.03 


0.31(0.30) 


0.14(0.12) 


0.09(0.07) 


0.11(0.09) 


0.34 


±0.03 


0.41(0.41) 


0.22(0.19) 


0.16(0.14) 


0.19(0.16) 


v2d y2 


0.18: 


t 0.025 


0.31(0.29) 


0.19(0.16) 


0.12(0.10) 


0.15(0.12) 


0.37 


±0.04 


0.40(0.39) 


0.27(0.25) 


0.20(0.17) 


0.23(0.20) 


vlh n/2 






0.08(0.06) 


0.06(0.05) 


0.05(0.04) 


0.04(0.03) 


0.25 


±0.07 


0.10(0.09) 


0.09(0.07) 


0.08(0.06) 


0.05(0.04) 


Orbit 








U8 Sn 












124 Sn 






Exp. 


PKA1 


PK03 


PKDD 


DD-ME2 


Exp. 


PKA1 


PK03 


PKDD 


DD-ME2 


v2d 5/2 


1.00 


±0.10 


0.91(0.93) 


0.86(0.89) 


0.88(0.90) 


0.92(0.94) 


1.00 


±0.09 


0.96(0.97) 


0.94(0.95) 


0.95(0.96) 


0.96(0.97) 


vlg 1/2 


0.78 


±0.19 


0.81(0.83) 


0.93(0.94) 


0.95(0.96) 


0.96(0.97) 


0.81 


±0.23 


0.93(0.94) 


0.96(0.97) 


0.97(0.98) 


0.98(0.98) 


v3s l/2 


0.80 


±0.10 


0.65(0.67) 


0.45(0.45) 


0.40(0.40) 


1.00(1.00) 


0.95 


±0.10 


0.88(0.90) 


0.77(0.78) 


0.73(0.74) 


1.00(1.00) 


v2d 3/2 


0.60 


±0.08 


0.61(0.62) 


0.52(0.52) 


0.46(0.46) 


0.56(0.58) 


0.75 


±0.10 


0.85(0.87) 


0.81(0.82) 


0.77(0.78) 


0.86(0.88) 


vl/in/2 


0.35 


±0.08 


0.17(0.16) 


0.20(0.18) 


0.20(0.18) 


0.13(0.10) 


0.38 


±0.12 


0.46(0.44) 


0.47(0.46) 


0.49(0.48) 


0.43(0.42) 



NL1 calculations the effective pairing interaction was taken as 
the Gogny force D1S with additional strength factor ( — 1.15 
whereas for the latter Gogny calculations it kept as the original 
one ATI . 



As seen from Table III the pairing energies are remarkably 
reduced in the calculations of PKA1, PK03 and PKDD with 
the optimized pairing force V^ f , which leads to weaker pair- 
ing contributions than those from NL1 and Gogny calcula- 
tions. In this work the pairing interaction is optimized as re- 
ferred to the OESs of binding energies of the odd isotopes. 
While in Ref. ED the strength factor f = 1.15 for NL1 is 
simply determined to reproduce the pairing energies of the 
Gogny calculations for Pb isotopes. Such deviations between 
the models may also originate from the inconsistence between 
the relativistic mean field and non-relativistic pairing poten- 
tial as mentioned before. As referred to the pairing energies, 
the PKA1 and PK03 calculations with original Gogny pairing 
force present slightly stronger pairing effects than the Gogny 
ones while PKDD, as well as DD-ME2, presents much weaker 
pairing effects. This discrepancy between the RHF and RMF 
models can be interpreted by the values of the non-relativistic - 
type effective mass [ 20 , 6 1 1 . With the inclusion of Fock terms, 
PKA1 and PK03 have a fairly large effective mass, close to 
the non-relativistic ones. Whereas in general the RMF mod- 
els present smaller effective masses, e.g., by PKDD and DD- 
ME2. Globally the level densities determined by PKA1 and 
PK03 are then larger than PKDD and DD-ME2 such that 
stronger pairing effects are obtained by the former ones with 
the same pairing interaction. This may also partially explain 
the reason why we have different pairing strength factor from 
the previous NL1 calculations. 

To provide a statistical understanding on the inconsistence 
between the non-relativistic Gogny pairing force and relativis- 
tic mean fields, Fig. |5]and Fig. [6]show the statistic behaviors 
of the deviations from the experimental data for both bind- 
ing energies and OESs, respectively. For comparison, the re- 




-5 5 10 15 20-5 5 10 

E - E (MeV) E - E (MeV) 



-5 5 10 15 



FIG. 5: (Color online) Distributions of deviations from experiment 
of the binding energies of 273 nuclei in the selected isotopes and 
isotones (1-MeV bins). Results obtained without (with) the effective 
pairing strength factor are shown in the upper (lower) panels. The 
best fitted Gaussian distributions [see Eq. |T0j] are also shown with 
black lines, and the corresponding parameters are listed in Table|IV| 




FIG. 6: (Color online) The same as Fig. |5]but for the OES [see Eq. 
{7}] of 213 nuclei with both even and odd nucleon numbers in the 
selected isotopes and isotones (0.1 -MeV bins). 
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TABLE III: The pairing energy of Pb isotopes obtained by the CDF calculations with PKA1, PK03 and PKDD. The results inside/outside 
the parenthesis are presented by the calculations taking the Gogny force D1S with/without pairing strength factor f as the effective pairing 
interaction. The corresponding results of RHB theory with NL1 with ( = 1.15 and of HFB theory with Gogny force |41 1 are also shown for 
comparison. 



A 


PKA1 


PK03 


PKDD 


NL1 




196 


-24 72('-17 76i 


-24 54f'-17 04 s ! 


- 1 7 47 (- 1 1 62) 






198 


-22.40(- 16.03) 


-22.40(- 15.38) 


-15.56(-10.14) 






200 


-19.40(-13.73) 


-19.52(-13.13) 


-13.14(-8.34) 






202 


-15.65(-10.83) 


-15.86(-10.28) 


-10.23(-6.26) 


-14.49 


-14.41 


204 


-11.14(-7.34) 


-11.34(-6.85) 


-6.85(-3.97) 


-10.41 


-10.49 


206 


-5.90(-3.50) 


-6.00(-3.11) 


-3.17(-1.48) 


-5.57 


-5.74 


208 


0.00(0.00) 


0.00(0.00) 


0.00(0.00) 


0.00 


0.00 


210 


-4.81(-2.93) 


-6.89(-4.30) 


-5.34(-3.66) 


-4.79 


-4.19 


212 


-8.58(-5.04) 


-12.00(-7.67) 


-9.61(-6.56) 


-8.66 


-8.64 


214 


-11.64(-6.37) 


-16.04(-10.33) 


-12.99(-8.87) 


-11.77 


-12.89 



TABLE IV: Parameters of the corresponding Gaussian distributions 
[see Eq. jT0}] for the binding energy in Fig. [5] and the OES in Fig. 
[6] Results obtained without/with the effective pairing strength factor 
are shown outside/inside the parenthesis. 



can be also demonstrated from the statistic variables shown in 





y>o 


x c 


ii' 


A 






Binding Energy 




PKA1 
PK03 
PKDD 


1.94(2.26) 
3.88(2.23) 
4.23(3.86) 


0.66(0.80) 
-0.50(0.03) 
-0.17(0.04) 


3.70(4.29) 
2.94(3.71) 
3.46(3.32) 


214.64(202.45) 
189.48(222.14) 
156.77(164.74) 








OES 




PKA1 
PK03 
PKDD 


0.89(1.57) 
0.95(1.54) 
1.48(1.70) 


0.19(-0.01) 
0.24( 0.03) 
0.08(-0.10) 


0.42(0.29) 
0.37(0.29) 
0.33(0.29) 


20.05(17.82) 
19.96(18.13) 
18.32(17.62) 



suits calculated with the original and optimized Gogny pair- 
ing forces are respectively shown in the upper and lower pan- 
els. In these two figures the solid curves denote the Guassian 
statistic fittings as, 



y = yo 



-2± 



(10) 



where yo is the minimum counts, A represents the area of the 
Gaussian distribution, and x c and w denote the statistic av- 
erages (ideally zero) and errors. In Fig. [5] and Fig. [6] it is 
clearly shown that the deviations of both binding energies and 
OESs from the data obey the Gaussian statistic behaviors. For 
the binding energy the deviations are not symmetrically dis- 
tributed on the negative and positive sides and more counts 
are found for the positive deviations, mainly due to the de- 
formation effects neglected in current calculations of Zr, Ce 
and Gd isotopes, as well as for the isotones. It is expected 
to get the results improved by including the deformation ef- 
fects, which in general lead to more bound groud states. For 
the OESs the statistic distributions are nearly symmetric, es- 
pecially for the PKA1 results with the optimized pairing force. 
Evident improvements due to the optimize pairing interaction 



Table IV It is seen that the statistic averages and errors on the 
OES deviations, i.e., the values of x c and w, are improved dis- 
tinctly. While the improvements on the binding energies are 
not so distinct and systematical as the OESs. Especially for 
PKA1 the statistic qualities even become worse, inconsistent 
with the evident improvements on the OESs. This may also 
indicate the inconsistence between the relativistic mean field 
and non-relativistic pairing interaction. 



IV. SUMMARY 

In summary, the neutron odd-even staggerings (OESs) of 
C, O, Ca, Ni, Zr, Sn, Ce, Gd and Pb isotopes, and the pro- 
ton ones of N = 50 and 82 isotones, are studied systemati- 
cally within the relativistic Hartree-Fock-Bogoliubov (RHFB) 
theory with PKA1 and PK03, and the relativistic Hartree- 
Bogoliubov (RHB) theory with PKDD and DD-ME2. By tak- 
ing the finite-range Gogny force D1S as an effective paring 
interaction, the neutron/proton OESs can be properly repro- 
duced by the selected effective Lagragians. 

Due to the inconsistence between the non-relativistic pair- 
ing force (Gogny D1S) and the relativistic effective La- 
grangians (PKA1, PK03, PKDD and DD-ME2), some sys- 
tematic discrepancies with the data of the OESs are found, 
namely, the OESs and coherently the pairing effects are some- 
what underestimated for light nuclei and overestimated for the 
heavy ones. Such inconsistence can be eliminated partially 
by introducing a Z- or A^-dependent strength factor into the 
pairing force, which present better agreements with the exper- 
imental data of the OESs, by improving the root-mean-square 
deviation cr rms of the OES about 10% for light nuclei and 30% 
for heavy ones. Similar improvements are also obtained on the 
description of single-particle configurations, especially in the 
calculations with PKA1. While with the optimization of pair- 
ing treatment, different even opposite systematical changes on 
the OESs and nuclear binding energies imply the remaining 
inconsistence between the relativistic mean fields and non- 
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relativistic pairing interaction, which may call for an unified 
relativistic mechanism in both p-p and p-h channels as future 
perspectives. 
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